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I. INTRODUCTION 



Owing to the non-perturbative nature of QCD in low energy regions, phenomenlogical 
models reflecting the characteristic of the strong interaction are widely used in the studying 
of the properties of hadrons, nuclear matter and quark matters [1-1 2]. They are based on 
different degrees of freedom, for example: nucleons and mesons, or quarks and gluons, in 
particular, on hybrid quarks and mesons. Some of these models have been proved to be 
successful. The quark mass density-dependent model(QMDD)[5] is one of such candidates. 

According to the QMDD model, the masses of u, d quarks and strange quarks (and the 
corresponding anti-quarks) are given by 

ruq^- — {i ^ u,d,u,d) (1) 

oUb 

B 

ms,s = msQ + - — (2) 

where is the baryon number density, m^o is the current mass of the strange quark, and 
B is the bag constant. At zero temperature 

nB = ^{riu + nd + ns), (3) 

where n„, n^, Ug represent the density of u quark, d quark, and s quark, respectively. The 
basic hypothesis Eqs.(l) and (2) corresponds to a quark confinement mechanism because if 
quark goes to infinite space, the volume of the system tends to infinite, ub approachs to 
zero and the rUq goes to infinite, and the infinite quark mass prevents the quark from going 
to infinite. The confinement mechanism is similar to that of the MIT bag model. 

Although the QMDD model can provide a description of confinement and explain many 
dynamical properties of strange quark matter, but it is still an ideal quark gas model and 
cannot explain the temperature T vs. density p deconfinement phase diagram of QCD and 
the properties of nuclear mattcr[13-14]. To overcome this difficulty, we have introduced a 
coupling between quark and nonlinear scalar field to improve the QMDD model in a previous 
paper [15]. We have found the wave functions of the ground state and the lowest one-particle 
excited states. By using these wave functions, we calculated many physical quantities such 
as root-mean-squared radius, the magnetic moment of nucleon to compare with experiments 
and come to a conclusion that this improved QMDD model is successful to explain the 



properties of nucleon. In ref[16], we extended this model to finite temperature and studied 
its soliton solution by means of the finite temperature quantum field theory. The critical 
temperature of quark deconfinement Tq and the temperature-dependent bag constant B{T) 
are found. The results of improved QMDD(IQMDD) model are qualitatively similar to that 
obtained from Frcidbcrg-Lee soliton bag model[17, 18]. 

Instead of studying the nucleon properties, we hope to employ the IQMDD to investigate 
the physical properties of nuclear matter in this paper. As was shown by the Walecka 
model[l] and the QMC model[7-10] early, a neutral vector field coupled to the conserved 
baryon current is very important for describing bulk properties of nuclear matter. The large 
neutral scalar and vector contributions have been observed empirically from NN scattering 
amplitude. The main qualitative features of the nucleon-nucleon interaction: a short range 
repulsion between baryons coming from a;-meson exchange, and a long-range attraction 
between baryons coming from cr-meson exchange must be included in a successful model. 
Obviously, if we hope to employ the IQMDD model to mimick this repulsive and attractive 
interactions, except the quark and cr-meson interaction, the uj meson and the qquj coupling 
must be added. This motivate us to introduce oj mesons and the qqoj coupling in the IQMDD 
model in this paper. In this new IQMDD model, the nonlinear scalar field coupling with 
quarks forms a soliton bag, and the qquj vector coupling gives the repulsion between quarks. 
We will prove that this model can give us a successful description of nuclear matter. 

The organization of this paper is as follows. In the next section, we give the main formulae 
of the IQMDD model under the mean field approximation at zero temperature. In the third 
section, some numerical results are contained. The last section contains a summary and 
discussions. 

II. FORMULAE OF THE IMPROVED QMDD MODEL 

The Lagrangian density of the IQMDD model is : 

L = V[^7'^^M -m, + g^a - giru;^]iP + ^d^ad^^a - U{a) - ^F^.F^'' + ^mluj^u;^ (4) 
where 

F^^ = d^ujy - dyUJ^ (5) 
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and the quark mass niq is given by Eqs.(l) and (2), rria and rriu) are the masses of a and cu 
mesons, and are the couphngs constant between quark-cr meson and quark-o; meson 
respectively. And 

U{c7) = ^m'y + ^-ba' + \ca' + B (6) 

-B= ^al + + '-at- (7) 

where is the absolute minimum of U{(t), U{(J^) = and U{0) = B. 

It can easily show the equation of motion for quark field in the whole space is 

Yl^{id, + gluj,) - (m, - gla)]^l: = (8) 

Under mean field approximation, the effective quark mass m* is given by: 

ml^mq- gla (9) 

In nuclear matter, three quarks constitute a bag, and the effective nucleon mass is obtained 
from the bag energy and reads: 

4 ^ r^l , AN 

M*N-^,E, = E,-nR'j^^l ^mf + k\^)dk (10) 

where quark degeneracy 7g=6, Kp is the Fermi energy of quarks. dNq/dk is the density of 
states for various quarks in a spherical cavity. It is given by [19]: 

N{k) = A{kRf + B{KRy + C{KR) (11) 



where 



A^'^. (12) 



^^m,^ - C(^) + (x)''''3.42(^ -15)2 + 100- ^^^^ 



C(^) ^ + (^ + ^)aTctan{—) - ^}. (15) 

^ A; ' 27r^3 ^ k V/ 2m/ ^ ' 
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Eqs. (12) and (13) are in good agreement with those given by multireflection theory[21, 22] 
and the Eqs. (14) and (15) are given by a best fit of numerical calculation for the MIT bag 
model. The curvature term C cannot be evaluated by this theory except for two limiting 

cases niq and rUq oo. Madsen[20] proposed the Eq. (15), but as was pointed out by 
Ref. [19], the best fit of numerical data is given by Eq. (14). This density of state has been 
employed by Refs. [13,14] to study the strangelets. 
The Fermi energy Kp of quarks is given by 

4 



where ub satisfies 



g-rfne (16) 



The bag radius R is determined by the equilibrium condition for the nucleon bag: 

= (18) 



5R 

In nuclear matter, the total energy density is given by 

2 1 11 

Sn^atter = // \l M*/ + pW + + 3^^^' + 3^^' + 4^^' (19) 

where 7Ar = 4 is degeneracy of nucleon, Kp is fermi energy of nucleon and is the density 
of nuclear matter 



In 



In Eqs. (19), g^j is the coupling constant between the nucleon and the cu meson and it 
satisfies — 3g^. As that of the QMC model[7], the a is yielded by the equation: 

2,2 S IN f^F Ml ,dMt^ 

mla + ha^ + ca^ = --i^ / d'p{—^)R (21) 



Eqs. (9)-(21) form a complete set of equations and we can solve them numerically. Our 
numerical results will be shown in the next section. 

III. NUMERICAL RESULT 

Before numerical calculation, let us consider the parameters in IQMDD model. As that 
of Ref.[l, 23], the masses of a;- meson and cr-meson are fixed as — 783 MeV, — 509 



MeV respectively. We choose the bag constant B — 174 MeV fm~^ to fit the mass of nucleon 
Mjv = 939 MeV. When B is determined, the parameters b and c are not independent because 
of Eq. (7). we choose the b is free parameter. There are still three parameters, namely, 
g^, g^, b are needed to be fixed in IQMDD model. 

To study the physical properties of nuclear matter, we investigate the nuclear saturation, 
the equation of state and the compressibility. The pressure of nuclear matter P is given by 



2 ^ ^matter 



P = Pb^ ^ (22) 

OPB PB 

where is the baryon density. The compressibility for nuclear matter reads: 

K = 9-^P (23) 
dpB 

at saturation point, the binding energy per particle E/A — —15 MeV, and the saturation 
density po — 0.15 fm~^. 

Our numerical results are shown in Fig. 1-Fig. 4. In Fig. 1, we choose ct;-meson and 
cr-meson satisfy a = 0, a) = and depict the bag energy as a function of bag radius at zero 
temperature. We find the stable radius of a "free" nucleon R — 0.85 fm. 

In Figs. 2-4 we show the effective mass M* of nucleon, the saturation curve and the 
equations of state of nuclear matter at zero temperature for IQMDD model respectively, 
where we fix the parameter b=-1460 (MeV), g^- — 4.67 and g^ — 2.44 respectively, and find 
E/A = -15 MeV and po = 0.15 fm"^ and K{po) = 210 MeV. We find our model can explain 
the properties of nuclear matter successfully. 

To illustrate our results more transparently, we show the dependence of the properties of 
nuclear matter on the parameters b,g^,g^ in Table. 1 for fixing binding energy E/A — —15 
MeV and po — 0.15 fm~^. We find that the compressibihty K{po) and effective nucleon mass 
Mj^(po) at saturation point all decrease when g^, g^ increase and b decreases. At was shown 
in Table. 1, the variational regions for K{pq) and M^(po) are small, and the decreasements 
of i^(po) and M^{pq) are slowly. 
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TABLE 1. Variation of the nuclear matter properties to b. 

b(MeV) gl gl K(po)(MeV) M^(po)(MeV) 



-800 4.59 2.35 218.8 782.9 

-1000 4.61 2.38 215.5 781.2 

-1200 4.64 2.40 213.6 778.5 

-1400 4.66 2.43 211.2 776.8 

-1600 4.69 2.46 208.1 774.3 

-1800 4.71 2.48 205.7 772.6 



It was pointed in the Refs. [24] early, adding a nonlinear scalar field in the model will 
cause unphysical behavior under mean field approximation in nuclear matter. This can 
easily be seen from Eq.(21) because the left hand side of Eq. (21) is a cubic order function 
of CT, and o" = is one of its solutions. There are two solutions in low-density regions. In 
Fig. 5 and Fig. 6, these two solutions are shown exphcitly for vs. ps curve and for M'^ 
vs. pb curve respectively where the parameters are fixing as 6 = —3655 (MeV). Noting that 
the term of non-linear scalar field is essential to form a soliton bag, we conclude that the 
unphysical branch cannot be avoided for the soliton solution under mean field approximation. 
Fortunately, the lower branch cannot be ended at the point (Mat = 939 MeV, pB — 0), and 
give us a experimental value of nucleon mass, we will give up this unphysical branch in our 
calculation. 

Finally, It is of interest to compare the properties of nuclear matter for IQMDD model 
and for the QMC model. Our results arc shown in Table. 2. we find their results are very 
similar. But as was pointed in our previous paper [15], the first advantage of the IQMDD 
model is that the MIT bag boundary constraint has been given up because it mimicks 
to a Friedberg-lee soliton bag model[13-16]. The second advantage of the IQMDD model 
is that the interactions between qqa and qquj are extended to the whole free space. We 
can easily write down the propagators of quarks, u-meson and a;-meson respectively and 
do the many-body calculations beyond mean field approximation. But for QMC model, 
the propagators of quarks, o"-meson and cj-meson cannot be written down easily because 
one must consider the multireflection by the MIT bag boundary as well as the effect of 
the interactions hmited into a nucleon space only. The IQMDD model provides a good 
substitute of the QMC model which is more suitable for the study of nuclear matter beyond 
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mean field. 



TABLE 2. Comparison of properties for IQMDD and QMC model. 







R(fm) 


9l 


9l 


K(po)(MeV) Mj^(po)(MeV) 


QMC 




0.80 


5.53 


1.26 


200 


851 


lQMDD(b= 


=0) 


0.85 


4.54 


2.21 


227 


798 


lQMDD(b= 


=-1460) 


0.85 


4.67 


2.44 


210 


775 


QHD-1 






Qa = 9.57 


= 11.6 


540 


522 



IV. SUMMARY AND DISCUSSION 

In summary, we present an Improved quark mass density dependent model which has the 
non- linear a meson field, and the uj meson field. The qqa coupling and the qquj coupling 
are introduced to mimick the repulsive and the attractive interactions between quarks in 
this model. It is shown that the present model is successful for describing the saturation 
properties, the equation of state and compressibility of nuclear matter. The effective nucleon 
mass decreases with baryon density in this model more rapidly than that of QMC model. 
After comparing the IQMDD model and the QMC model, we come to a conclusion that the 
IQMDD model is a good substitute for QMC model. 
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FIG. 1: The bag energy as a function of bag radius at zero temperature for w = 0, w = 0. 
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FIG. 2: Effective nucleon mass vs. baryon density at zero temperature where tlie parameters 
gl = 4.67, = 2.44,6 = -1460 (MeV). 
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FIG. 4: Pressure of nuclear matter as a function of pB- the parameters is same as that of Fig. 2. 



13 




70- 

t ■ 

't> 60- 

50- 




0.00 0.05 0.10 0.15 

Fig.5 



FIG. 5: the a field vs. baryon density for b=-3655 (MeV), gl = 5.23, = 3.12. 
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FIG. 6: Effective nucleoli mass M* vs. baryon density, the parameters is same as that of Fig. 5. 
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